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Introduction 

Very  few  problems  involving  scattering  of  a  wave  by  an  obstacle  can  be 
solved  exactly.  The  object  of  the  present  paper  is  to  examine  a  method  for 
numerical  solution  of  the  integral  equation  for  diffraction  by  an  obstacle 
for  which  solutions  are  available,  namely  a  strip,  with  a  view  to  the  even¬ 
tual  application  of  the  method  to  scattering  by  obstacles  of  more  general 
shape. 

We  consider  the  problem  of  scalar  two-dimensional  diffraction  of  the 
incident  wave 


ik(x  cos  0+  y  sin  ©) 

-  e  , 

by  a  soft  strip  lying  in  y  =  0,  -a  S  x  §  a.  It  is  requirelto  solve  the 
steady-state  wave  equation 

IF  +  ~W  +  k^t  =  0  ' 

where  <t>,  =  <)>.  +  $,  is  the  total  velocity  potential  consisting  of  the  sum  of 
1  1 

the  incident  acid  scattered  potentials  <K  and  <t>  respectively.  The  scattered 
potential  must  satisfy  a  radiation  condition  at  infinity.  The  boundary  con¬ 
dition  on  the  strip  is 


or 


=  °* 

.  ik  x  cos  © 

<t>  =  -e 


y=0,  -aSxSa, 
y  =  0,  -a  S  x  S  a  . 


(1) 


By  superposition  of  simple  sources  of  strength  f(x)  per  unit  length  on  the 
strip,  the  scattered  potential  at  any  point  (x,y)  of  space  is  given  by 
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If  we  let  (x,y)  tend  to  a  point  on  the  strip  we  obtain,  on  applying  the 
boundary  condition  (l),  an  integral  equation  for  the  unknown  function  f(x): 


-  |  f  f(6)  H^1  ^  (k  lx  -  ||)  d£  =  elkx  cos  0,  (-a  *  x  S  a).  (3> 

“  -a 

Of  the  voluminous  literature  on  this  subject  we  quote  the  exact  Mathieu 
function  solution  of  Morse  and  Rubinstein  (7),  the  low-frequency  expansions 
of  Millar  (5)  using  an  integral  equation  method  of  Bouwkamp,  and  the  high- 
frequency  asymptotic  expansions  of  Millar  (4). 

Low-frequency  expansions  have  been  obtained  by  a  variational  technique, 
using  analytical  methods,  by  De  Hopp  (3).  An  estimate  of  the  scattering 
cross-section  for  normal  incidence,  using  a  crude  constant  approximation  for 
the  source-strength  f(x)  on  the  strip,  has  been  given  by  Erdelyi  and  Papas  (2). 
The  present  paper  can  be  regarded  as  an  extension  of  the  method  of  Erdelyi 
and  Papas  using  the  method  described  in  (8).  The  difficulty  in  applying  a 
variational  principle  to  scattering  by  obstacles  with  sharp  edges  is  that 
the  unknown  functions  have  singularities  at  the  sharp  edges.  If  we  attempt 
to  use  trial  functions  with  the  correct  singularities  at  the  edges  we  find 
that  it  is  difficult  to  perform  the  necessary  integrations.  The  integra¬ 
tions  can  be  performed  for  constant  trial  functions,  but  these  do  not 
simulate  the  required  singularities.  In  this  paper  we  adopt  the  following 
compromise.  To  simulate  the  function 

f(x)  =  1  /\/ ( a2  -  x2), 


we  use  the  approximation 


f(x)  =  ja  -p  <  x  <  p  , 

Ip  p  <  |x|  <  a  , 


3 


where  a,  P,  p  are  determined  from  the  variational  principle.  It  is  found  by 

the  method  described  later  that  typical  values  for  normal  incidence  and  low 

frequency  are  0/a  =  p  =  0.92. 

The  variational  principle 

Consider  the  variational  expression 
pP  pP 

l(E,F)  =  /  E(x)  g(x)  dx  +  /  F(x)  h(x)  dx 

^  a  '-'a. 

pP  pP 

-  /  /  K(x,|)  E(x)  F(|)  dx  d|  ,  (4) 

J  a.  v  a 

where  g,  h,  K  are  known  functions,  and  E,  F  are  to  be  determined  so  as  to 
extremise  l(E,F).  We  set 

E(x)  =  e(x)  +  St|(x),  F(x)  =  f(x)  +  e£(x)  ,  (5) 

where  we  assume  that  e,  f  are  the  functions  which  we  wish  to  find,  namely  the 
functions  which  extremise  l(.E,F).  In  the  usual  way  T)(x)>  £(x)  are  arbitrary 
functions.  On  varying  the  parameters  8,e  the  functions  E,  F  vary  round  the 
required  functions  e,  f.  On  inserting  (5)  in  (4)  we  see  that 

l(E,F)  =  l(e,f )  +  A  +  B, 

where  A  is  first  order  in  5  and  e,  and  B  is  second  order: 
pP  pP 

A  =  e  /  £(x)  {h(x)  -  /  K(|,x)  e(|)  d|)  dx 

^  a  J  a 

+  5  r  fi(x)  {g(x)  -  f  K(x,g)  f(g)  d£}  dx. 

'-'a  J  a 


If  I  is  stationary  for  variations  pf  e,  5  round  zero,  the  coefficients  of 
e,  8  are  separately  zero.  Since  £(x),  t](x)  are  arbitrary,  e(x),  f(x) 

satisfy  the  integral  equations 
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K(x,g)  f(g)  dg  =  g(x)  , 


(a  i  x  g  p)  , 


K(§,x>  e(g)  dg  =  h(x)  , 


(a  §  x  g  p) 


If  e,  f  must  satisfy  these  equations  then  we  see  from  (4)  that 

pP  pP 

l(e,f)  =  /  e(x)  g(x)  dx  =  /  h(x)  f(x)  dx  . 


Suppose  that  we  can  guess  the  shape  of  e(x),  f(x)  so  that  we  can 
approximate  e(x),  f(x)  hy  the  expressions 

E(x)  =  c'~r(x)  ,  F(x)=dV(x), 

where  £-  /d'  are  known  functions,  and  C,  D  are  constants  which  have  to  he 
determined.  If  these  expressions  are  substituted  in  (4)  we  can  determine 

aj  ^  j 

optimum  values  for  C  and  D  hy  setting  This  gives 


q  _ _ p£-;_hJ'J 


aTsTV) 


I(E,F)  - 


a  (a >,#r  ’ 


where  we  use  the  notation 


(p,q.)  =  /  P(x)  q(x)  dx, 


A(p,q)  = 


K(x, g )  p(x )  q(g )  dx  dg. 


a  w  a 


Application  of  the  variational  principle  to  diffraction  hy  a  strip 

We  now  apply  the  results  obtained  in  the  last  section  to  the  problem  of 
diffraction  by  a  strip.  By  using  the  asymptotic  expansion  of  th6  Hankel 
function  in  (2)  it  is  readily  seen  that  the  far-field  scattered  by  the  strip 


is  given  by 
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A(e,e)  (&&r)-4  .l(kr+l”)  , 


where 


A  (0,0)  =  - 


•(s> 


a-ikgcos0 


d|. 


(9) 


-a 


In  these  expressions  we  use  polar  coordinates  (r ,0)  and  write  A (0,©)  to 
remind  us  that  A  depends  on  both  the  angle  0  at  which  we  observe  the  diffracted 
wave,  and  the  angle  0  of  incidence  of  the  incident  wave  (since  f(g)  depends  on 
©).  As  shown  by  several  authors  (for  instance,  de  Hoop  (5))  the  scattering 
cross-section  per  unit  width  of  the  strip  for  a  wave  incident  at  angle  ©  is 


given  by 

t(0)  =  -  (l/4ka)  Im  (A(0,©))  .  (10) 

We  obtain  a  variational  principle  for  A (0,©)  by  first  of  all  identifying 

the  integral  equation  (5)  with  (6a).  This  indicates  that  f(|)  can  be  taken  as 

the  same  function  in  both  equations  and  we  can  set 

K(x,£ )  =  -  (kix  -  ||),  g(x)  =  elkx  cos  0.  (ife) 

If  now  we  identify  minus  the  expression  in  (9)  with  (j)  we  see  that 

h(x)  =  e"ikX  COS  0  ,  (HD) 

and  from  (6b)  the  function  e(x)  must  satisfy  the  integral  equation 

-  Ji  «(S>  h5,’ >  <k  I  ^  -  ||)  a5  .  e-ito:  ».  (12) 

-a 


From  (7),  (9)  and  the  above  identifications  we  see  that  the  exact  value  of 
A(e,0)  is  equal  to  -l(e,f)  so  that  from  (8)  the  variational  expression  for 
A (0,©)  is 


A(6>,0) 


ikx  cos  © 
e 


dx 


a 


6^(x) 


-ikx  cos 
e 


dx 


£  \)t>(x)^(i)  dx  d| 
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This  is  the  form  of  the  variational  principle  used  by  other  writers  on  this 

subject,  for  example  Morse  and  Feshbach  (6)  and  de  Hoop  (3). 

For  our  purposes  it  is  more  convenient  to  proceed  in  the  following  way. 

Consider  the  integral  equations 
a 

K(l  x  -  |  J.)  f(g)  a£  =  g(x)  ,  (-a  g  x  g  a),  (13a) 

•a 


7 


In  an  exactly  similar  way,  equation  (13!)  can  be  replaced  by  the  following 
pair,  using  notations  analogous  to  (l4): 

f  K  (x,|)  e  (g)  dg  =  h  (x)  ,  (17a) 


K_(x,g)  e_(| )  d|  =  h_(x). 


(17b) 


On  comparing  (l6a),  (17a)  with  (6a,b)  we  obtain  from  (4)  a  variational 
expression  which  we  shall  denote  by  I+  (E+,  F+)  where  E+  and  F+  are  approxi¬ 
mations  to  e  and  f  .  The  exact  value  of  this  expression,  from  (7),  (l4),  isT 


I+(e+,f+)  =  /  e+(x)  g+(x)  dx  =  /  h+(x)  f+(x)  dx 


=  4  /  (e(x)  +  e("x))  (g(x)  +  g("x))  ^ 


/a 

(e(x)  g(x)  +  e(x)  g( -x ) )  dx. 


(l8a) 


Similarly  from  (l6b),  (17b)  we  can  derive  a  variation  expression  which  we 
denote  by  I  (E__,F  )  where  E_,  F_  are  approximations  to  e_,  f_.  The  exact  value 
of  this  expression  is  readily  shown  to  be 


r\  a 

I_(e_,f_)  =  l  (e(x)  g(x)  -  e(x)  g(-x)}  dx 


(18b) 


On  adding  (l8a,b),  and  remembering  that  A (0,0)  is  equal  to  -l(e,f)  whose 

value  is  given  in  (j )  we  see  that 

A(e,e)  =  -2  (I+(e+,f+)  +  I_(e_,f_ ))  .  ( 

The  variational  expression  for  I+  is 
n  a  n  a 

Ij.  =  /  E+(x)  g+(x)  dx  +  /  F+(x)  h+(x)  dx 


-  J  J  K+(x,|)  E+(x)  F+(|)  dx  dg  , 


(20a) 
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where  from  (l"l) 


g, (x)  =  cos  (kx  cos  ©),  h,(x)  =  cos  (kx  cos  0). 

T  4' 


(20b) 


From  this  point  onwards  it  is  convenient  to  set  a.  =  1 .  We  also  substi¬ 
tute  k  for  k  to  remind  us  of  this  assumption. 


We  take  the  trial  functions  F^j  E+to  be  st ep- fun ct ions : 

F+(x)  =  |  a,  -p  <  x  <  p  E+(x)  =  j'  7,  -q  <  x  <  q  , 


j  P  <  (xl  <  1 

Substitution  in  (20a)  gives 


!  5>  q  <  ix  I  <  1  . 
1 


I+  =  7D1  +  SD2  +  aCi  +  (3C2  -  ttyA-n  ~  Qi8A12  -  P/A21  -  |35A22  (22) 


where  from  (20b) 


t  s  ,  sin  cos 

g,  (x)  dx  =  - - - 

+  k  cos  0 


D~  =  /  a  (x)  dx  =  sin  cos0)  -  sin  (k  Q  cos  ®) 
Jq  +  £  cos  0 


Similarly 


sin  (sip  cos  0, 
K.  cos  0 


>  ^2 


(k.  cos  0)  -  sin  (r.p  cos  0 ) 
it.  cos  0 


pP  p 1  pP  4 

An  =  /  /  K  (x,g)  dx  dg  ,  A12  =  /  / 

v  Q  d  o  d  o  d  q 


K+(x,| )  dx  d£  , 


pi  n  q  pi  pi 

Agi  =  /  /  K+(x,|)  dx  dj  ,  A22  =  /  / 

OpJ0  p  ^  q 


K+(x,|)  dx  d|  . 


When  K+(x,|)  has  the  special  form  given  in  (l6c)  these  double  integrals 
can  be  simplified  by  means  of  the  following  results  (cf.  Noble  (8)): 

n  Td  n  d 

j  J '  K(  I  x  -  |  |)  dx  d|  =  J((b  -  c  I)  +  J(|a  -  d  |)  -  j(|b  -  d  |)  -  J(  |a  -  c  |), 


b  nd 


K(x  +  g)  dx  dg  =  J(b  +  d)  +  j(a  +  c)  -  J(b  +  c)  -  j(a  +  d). 


9 


(23) 
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where  A  =  AX1  A22  -  A12  A21. 

Next  consider  the  variational  expression  I  derived  from  (l6b),  (l7h): 
n  a  p  a 

I  =  /  E  (x)  g  (x)  dx  +  /  F  (x)  h  (x)  dx 


-  /  /  K_(x,|)  E_(x)  F_(g)  dx  d£  ,  (24) 

O  o 

where 

g_(x)  =  i  sin  (kx  cos  0)  ,  h  (x)  =  -  i  sin  (kx  cos  0). 

We  again  set  a  =  1 ,  replace  k  "by  and  take  the  trial  functions  E_  to  he 

odd  step-functions : 


-Oi, 

-1  <  X  <  -p, 

e_(x)  =  ;-p, 

-1  <  X  < 

-q  • 

'  0, 

I 

pq  <  v> 

°> 

|X|< 

q  > 

(25) 

i  a. 

p  <  x  <  1 

1  p. 

q  <  x  < 

i  . 

Then  (24)  becomes 


where 


I_  =  PP  +  aQ  -  c*pB 

_  i  fcos  ( kp  cos  0)  -  cos yc  cos  0) 


Q  =  - 


K  COS  >0 

i  f cos  (.fcp  cos  &)  -  cos  (k  cos  6 ] 
'k  cos  0 


1  1 

B  =  f  f  K  (x,  £ )  dx  d£ . 

J  p  J  q 

By  following  the  procedure  explained  previously  in  connection  with  the  A. . 

i  J 

we  can  show  that 

41k2B  =  L(b(1  +  p»  +  L(k(l  -  p)}  +  L(k(  1  +  q)}  +  L{k(1  -  q)) 

-  L{n(p  +  q)}  -  L{klp  -  q |)  -  L(2.k). 
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On  setting  the  derivatives  of  I  with  respect  to  a,  (3  equal  to  zero  we  find 
the  optimum  values  of  these  parameters,  and  substitution  in  (2 6)  gives 

I_  =  pP  =  aQ  =  PQ/B.  (27) 

4  Numerical  results  for  the  scattering  cross-section 

The  numerical  results  reported  below  are  confined  to  the  scattering 
cross-section  per  unit  width  of  the  strip,  t(©),  defined  in  (10)  above. 

We  then  require  only  A(©,0)  so  that  we  can  set  9  =  0,  P  =  q,  in  the 

expressions  in  the  last  section.  On  using  (lO)  and  (19);  remembering  that 
we  now  set  a  =  1,  k  =  K.  in  (10),  we  have 

t(@)  =  ~  Im  (I  +  I  )  =  Re(j  +  J  ),  say, 

C.K  T  -  T 

whece  we  can  show  from  (23  )>  (27)  that 

j+  =  -  i+ 


where  we  have  separated  out  the  dominant  part  of  J  ,  and 

T 

M(it,p)  =  L{k(1  -  p)}  -  L(k(1  +  p)}.  Also 
J_  =  -  i(i/'K)  1  _ 

_ _ 2fcos  (  KP  cos  ©)  -  cos  (*K  cos  &)~\ 

~  k  cos^Q  (2[L{  n(1+p)}  +  L{k(1  -  p)]]  -  L(2k)  -  L(2icp)] 

The  optimum  value  of  p,  which  is  so  far  undetermined,  is  found 
empirically  by  using  the  variational  property  of  the  two  above  expressions. 
Thus  for  a  given  value  of  k.  and  0,  the  quantity  J,  is  evaluated  for  various 

-r 

p,  and  the  required  value  of  p  is  that  for  which  J+  is  stationary.  A  diffi¬ 
culty  arises  since  J+  is  complex  and  the  optimum  value  of  p  may  not  be  the 
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same  for  the  real  and  imaginary  parts.  Fortunately  it  turns  out  that  the 
optimum  values  of  p  are  not  essentially  different  when  k  is  not  too  large., 
as  illustrated  in  a  numerical  example  below.  The  optimum  value  of  p  need 
not  of  course  be  the  same  for  J+  and  J  since  these  are  obtained  from  inde¬ 
pendent  integral  equations. 

The  above  calculations  to  determine  the  optimum  value  of  p  were  carried 
out  for  various  angles  of  incidence  ©  =  0(10)90  degrees  and  k  =  0.4(0.4)2.8. 
The  results  were  similar  for  all  angles  of  incidence  and  figures  for  a 

typical  example,  namely  ©  =  50° >  are  given  in  Tables  1-4,  For  small  k.  the 

optimum  value  of  p  for  both  the  real  and  imaginary  parts  of  J+  was  approxi¬ 
mately  0.92  as  shown  in  Tables  1  and  3.  The  optimum  value  increases 
slightly  as  k  increases  to  2.0,  and  then  increases  fairly  rapidly  to  unity 
as  k  increases  from  2.0  to  3-0-  The  optimum  value  of  p  for  the  imaginary 
part  of  J+  increases  slightly  more  rapidly  than  the  optimum  value  for  the 
real  part  but  no  significant  error  is  involved  in  the  assumption  that  p  is 

a  constant  equal  to  0.92  for  all  k  between  0  and  2.0.  For  small  k  the 

optimum  p  for  both  the  real  and  imaginary  parts  of  J_  was  about  0.68.  The 
optimum  value  for  the  imaginary  part  increases,  and  for  the  real  part  it 
decreases,  as  p  increases,  but  the  estimates  of  the  real  and  imaginary 
parts  of  J  are  not  seriously  in  error  if  we  assume  p  constant  equal  to 
0.68  for  all  k  between  0  and  2.0. 

Some  typical  results  for  t(©)  are  plotted  in  Figures  1  and  2,  where  the 
optimum  values  of  p  for  the  real  parts  of  J+  and  J_  have  been  used.  In 
Figure  1  our  results  for  t(©)  are  plotted  against  k  for  k  =  0.4  to  2.0,  and 
0  =  0,  30,.  60°.  Also  in  Figure  1,  for  the  sake  of  comparison,  are  the  exact 
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TABLE  1 


Values  of  Re[J+]  for  0  =  50%  k'  =  0.4  to  2.8,  p  =  0. 


to  O.96 


■ 

0.88 

. 

1 

0.90 

0.92 

0.94 

0.96 

0.4 

1  1.1278 

1.1284 

mm 

0.8 

0.8824 

0.8829 

SB 

1.2 

0.7362 

0. 73.65 

0.7567 

0.7367 

0.7362 

1 .6 

0.6094 

0. 6095 

0.6097 

0.6098 

0.6097 

2.0 

0.4933 

0. 4932 

0. 4930 

1  0. 4930 

0. 4930 

2.4 

0.3932 

0.3924 

0.3916 

0.3907 

0.3898 

2.8 

0.3161 

0.3143 

0.3121 

0.3094 

0.3060 

TABLE  2 

Values  of  Re[J  ]  for  0  =  50°,  k  =  0. 4  to  2.8,  p  =  0.62  to  0.70 


K\ 

0.62 

0.64 

Q. 66 

0. 68 

0.70 

0.4 

0. 00559 

0.00561 

0. 00563 

0.8 

0.0350 

O.0351 

0. 0352 

0.0352 

0.0351 

1 .2 

0. 0909 

0.091 1 

0.0912 

0.0912 

0.0910 

1.6 

0. 1 656 

0.1659 

0. 1659 

0.1657 

0.1652 

2.0 

0.2486 

0.2488 

0.2485 

0.2479 

0.2468 

2.4 

0.3292 

0.3289 

0.3279 

0.3262 

0.3238 

2.8 

0.3970 

0.3950 

0.3920 

0.3877 

0.3820 
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TABLE  3 


Values  of  Im[J  ]  for  0  =  50° ,  k  =  0.1!  to  2,0,  p  =  0.88  to  0. 96 


B 

0.88 

0.90 

0.92 

0.94 

0.96 

0.4 

-1.2512 

-1.2513 

-1.2513 

-1 .2512 

-1 .2509 

0.8 

-0.5909 

-0.5906 

-0.5905 

-0.5905 

-0.5907 

1 .2 : 

-0.3078 

-0.3074 

-0.3071 

-0.3070 

-0.3073 

1.6 

-0.1541 

-0.1536 

-0.1532 

-0.1529 

-0. 1529 

2.0 

-0.0741 

-0.0737 

-0. 0733 

-0. 0728 

-0. 0724 

TABLE  4 

Values  of  Im[Jj  for  0  =  50  ,  k  =  0.4  to  2.0,  p  =  0.62  to  0.70 


\ 

\ 

\P 

K  \ 

0.62 

_  _ 

0.64 

0.66 

0. 68 

0.70 

0.4 

-0.1080 

-0. 1 082 

-0.1083 

-0.1084 

-0. 1 084 

0.8 

-0. 1918 

-0. 1 922 

-0.1924 

-0.1926 

-0. 1 925 

1 .2 

-0.2486 

-0.2492 

-0.2496 

-0.2499 

-0.2500 

1.6 

-0.2779 

-0.2788 

-0.2796 

-0.2803 

-0.2808 

2.0 

-0.281  4 

-0.2830 

-0.2846 

-0.2862 

-0.2878 
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solutions  of  Morse  and  Rubinstein  (t)  for  0  =  0%  and  the  curves  of  Millar 
(4)  and  (5).  It  should  be  noted  that  Millar's  curves  consist  of  a  low- 
frequency  approximation  extending  up  to  k  ~  1 ,  and  a  high-frequency  approxi¬ 
mation  extending  down  to  is  £  1.25#  These  two  approximations  do  not  match 
exactly.  In  Figure  2  our  results  for  t(©)  are  plotted  against  k  going  from 
k  =  0.4  to  4.0  for  normal  incidence  (©  =  90° )•  Also  included  are  the  exact 
solutions  of  Morse  and  Rubinstein,  Millar's  solutions  for  normal  incidence, 
and  the  solution  of  Erdelyi  and  Papas  (2).  Our  results  agree  well  with  the 
exact  solution  of  Morse  and  Rubinstein  for  0. 4  <  k  <  5  in  Figure  2  and 
0. 4  <  k  <  2  in  Figure  1 .  Our  results  are  also  in  fair  agreement  with 
those  of  Millar,  the  agreement  being  worst  at  grazing  incidence  where  in 
any  case  Millar's  low-  and  high-frequency  approximations  do  not  join 
smoothly  in  the  region  1  <  k  <  1.25. 

The  solution  of  Erdelyi  and  Papas  for  normal  incidence  given  in 
Figure  2  is  obtained  from  the  variational  principle  which  is  used  in  this 
paper  but  a  trial  function  is  assumed  which  is  simply  constant  over  the 
whole  strip.  The  trial  function  of  Erdelyi  and  Papas  is  obtained  by 
taking 

p  =  q  =  1,  F+(x)  =  E+(x)  =  a,  (-1  <  x  <  1 ) 

in  equation  (21 )  above.  We  have  already  mentioned  that  the  optimum  value 
of  p  for  the  even  solution  tended  to  unity  as  k  increases  from  2.0  to  5.0, 
so  that  the  graph  of  our  results  in  Figure  2  tends  to  the  solution  of 
Erdelyi  and  Papas  for  is  >  5.0.  Unfortunately  the  asymptotic  behaviour 
of  the  Erdelyi  and  Papas  solution  does  not  agree  with  that  of  the  exact 
solution  as  k  tends  to  infinity.  However  our  results  show  a  considerable 


Sceitt&riiuj  crass^ed\On  i  (©)  0.S  <x  -f u/^ction  o-\ 
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improvement  over  those  of  Erdelyi  and  Papas  for  k  less  than  about  3.0. 

The  trial  functions  used  in  this  paper  are  quite  crude  in  that  they 
are  step -functions  with,  effectively,  only  one  step.  In  spite  of  this 
the  results  are  remarkably  good.  It  is  also  noteworthy  that  empirically 
it  appears  that  no  great  loss  of  accuracy  is  involved  in  computing  the 
scattering  cross-section  if  it  is  assumed  that  for  all  0  we  have 
p  =  q  =  0.92  for  the  even  function  (21 )  and  p  =  q  =  0.68  for  the  odd 
function  (25),  in  the  range  0  <  k  <  2.5.  It  is  the  chief  conclusion  of 
this  paper  that  good  results  can  be  obtained  with  these  simple  trial 
functions.  This  encourages  the  hope  that  it  may  be  possible  to  apply 
the  same  method  to  calculate  the  scattering  by  objects  of  more  compli¬ 
cated  shapes  in  the  awkward  wavelength  region  where  the  wavelength  is 
comparable  with  the  size  of  the  obstacle. 
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